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Abstract: We study analytically quasinormal modes in a wide variety of black hole spacetimes, including 
d-dimensional asymptotically flat spacetimes and non-asymptotically flat spacetimes (particular attention 
has been paid to the four dimensional case). We extend the analytical calculation to include first-order 
corrections to analytical expressions for quasinormal mode frequencies by making use of a monodromy 
technique. All possible type perturbations are included in this paper. The calculation performed in this 
paper show that systematic expansions for uncharged black holes include different corrections with the 
^ . ones for charged black holes. This difference makes them have a different n-dependence relation in the 
first-order correction formulae. The method applied above in calculating the first-order corrections of 
quasinormal mode frequencies seems to be unavailable for black holes with small charge. This result 
supports the Neitzke's prediction. On what concerns quantum gravity we confirm the view that the In 3 
in d = 4 Schwarzschild seems to be nothing but some numerical coincidences. 
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1. Introduction 



The stability of black holes has been discussed since researchers found that black holes were shown to 
radiate and evaporate when we add in the ideas of quantum mechanics to them|l], |2| and hence people 
tried to make sure whether the black hole solution under consideration was really a stable one of the 
classical theory. The pioneering work on this problem was carried by Regge and Wheeler Q, who focused 
on analyzing the linear stability of four dimensional Schwarzschild black hole. They found that one can use 
a Schrodinger-like equation to describe the linear perturbations. This work was latter extended to many 
other black hole solutions and is now known as the quasinormal modes (QNMs) which can be described as 
a "characteristic sound" of black holes ( A lot of investigation have been made on this subject ||, |{| |7|). 
QNMs are excited by the external perturbations(may be induced, for example, by the falling matter). 
They appear as damped oscillations described by the complex characteristic frequencies which are entirely 
fixed by the parameters of the given black hole spacetime, and independent of the initial perturbation ||. 
These frequencies can be detected by observing the gravitational wave signal || : this makes QNMs be of 
particular relevance in gravitational wave astronomy. 

Although the QNMs are important in the observational aspects of gravitational waves phenomena 
mentioned above, there is suggestion that the asymptotic QNMs may find a very important place in Loop 
Quantum Gravity (LQG). Recently Hod [|l]] made an interesting proposal to infer quantum properties of 
black holes from their classical oscillation spectrum . The idea was based on the Bekenstein's conjecture 
[11] that in a quantum theory of gravity the surface area of a non-extremal black hole should have a discrete 
eigenvalue spectrum. The eigenvalues of this spectrum are likely to be uniformly spaced. According to the 
numerical values computed by Nollert Q and later confirmed by Andersson [12|, Hod observed that the 
real parts of the asymptotic form of high overtones of a Schwarzschild black hole can be written as : 

= (2n+l)vri + ln3, 

where Tjj is the Hawking temperature. Using this conjecture and Bekenstein's conjecture, he obtained 
the Bekenstein-Hawking entropy for the Schwarzschild black hole. He also showed that this approach 
is compatible with the statistical mechanical interpretation of black hole entropy. Later, Dreyer |l3|] 
noted that this is helpful for the calculation of the so called Barbero-Immirzi parameter, a free parameter 
introduced as the Barbero-Immirzi connection in the calculation of LQG. The only way so far that could 
be used to fix this parameter comes from black hole entropy. Dreyer use the conjecture made by Hod, he 
fixed the value for the Immirzi parameter by ln ^- . Using this value, he suggested that "the appropriate 
gauge group of quantum gravity is SO (3) and not its covering group SU (2)". However, more recently 
Corichi[14] argued that the LQG allows us to keep SU (2) as the gauge group, and at the same time have 
a consistent description with the results of Dreyer. He reconsidered the physical process that would give 
rise to the quasinormal frequency (QN frequency) as mentioned in [13]: an appearance or disappearance 
of a puncture with spin jrnin- Taking into account of the local conservation of fermion number, Corichi 
obtained that "the minimum allowed value for the 'spin' of the resulting free edge is j m in = 1" : this agrees 
with the results of Dreyer. However, all this problems are far from being resolved. 

As mentioned in the above paragraph, the Hod's conjecture was based on the numerical results evalu- 
ated by Nollertpf]. People may think this agreement is just a coincidence. Motl[j^, M, however, confirmed 
analytically Nollert 's result by two different methods. In Ref.[15], the author used Nollert 's continued 
fraction expansion for the 4-dimensional Schwarzschild and showed that the asymptotic QN frequencies 
are in good agreement with Hod's result. The monodromy technique was first introduced in [|l(| to an- 
alytically compute asymptotic QN frequencies and later extended in so that it can also be used in 
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the computation of d-dimensional asymptotically dS, AdS spacetime. A question, however, produced as 
showed in [16, 17] is that the suggestion in [p], O] was proved not to be universal and be only applicable 



to the Schwarzschild solution. It is necessary to stress that even if the ideas in 1C, 13] turn out not to be 
universal, it is still the case that QN frequencies will play some role in the realm of quantum gravity, since 
the LQG is far from being successful, and recent studies show that they have interpretation in conformal 
field theory through AdS/CFT correspondence [18, IS, 20, 21]. Moreover, The possible appearance of ln2 



in the asymptotic frequencies |16| could support the claims [Q that the gauge group of LQG should be 
SU(2) despite the In 3 for Schwarzschild. 

It is important to note that the behavior of QNMs for the lowest modes (frequencies with a smaller 
imaginary part) is totally different from the one for the high overtone. We have shown that the real part 
of the QN frequencies approaches some certain value as the imaginary part approaches infinity. There has 
no such phenomena for the low- lying modes. A question of particular relevance that immediately follows 
is what will happen for the QNMs when the imaginary part has a middle value between these two extreme 
cases, and why these two cases behavior so differently (the resolution of this problem is helpful to deduce 
analytically the asymptotic value of the QN frequencies). Recently Musiri and Siopsis 



22], ^3[ have studied 

in detail about this question for Schwarzschild in asymptotically flat and asymptotically AdS spacetime. 
Extending the technique introduced in []hJ] to obtain a systematic expansion including corrections in 
1/y/u), they obtained the j and I dependence of the first correction for arbitrary j. Their results are in 
good agreement with the results obtained by numerical methods in the case of scalar and gravitational 
waves. However, they have discussed nothing about more general spacetime background, such as the case 
of Reissner-Nordstrom (RN) black holes and higher dimensions. The main purpose of this paper is to study 
analytically the first-order correction to the asymptotic form of QN frequencies for more general sapcetime 
background. In this work we shall make use of the remarkable results obtained by Ishibashi and Kodama 
(we refer the reader to j 



2q , [26! f° r detail on this subject). They studied in detail the perturbation 



theory of static, spherically symmetric black holes in any space-time dimension d > 3 and allowing for the 
possibilities of both electromagnetic charge and a background cosmological constant. According to them, 
the perturbations come in three types: tensor type perturbations, vector type perturbations and scalar 
type perturbations, and linear perturbations in d-dimensions can be described by a set of equations which 
may be denoted as the Ishibashi-Kodama (IK) master equations. 

The organization of this paper is as follows. In next section we apply the monodromy method in- 



troduced in [16] to the IK master equations expanded near the several singularities in the complex plane 
to analytically compute both zeroth-order and first-order asymptotic quasinormal frequencies for static, 
spherically symmetric black hole spacetimes in dimension d > 3. Section 3 is the last section of the paper, 
where we have a discussion about our results, listing some problems encountered in this paper. Some future 
directions are also included in this section. The last section is our appendices. In appendix A we make 
use of expanding the tortoise coordinate to first order at the singularities in the spacetimes considered in 
this paper, providing a full analysis of the potentials at several singularities in the complex plane, and 
obtaining a list of first-order IK master equation potentials. 



2. Perturbative Calculation of Quasinormal Modes 

In this section we first review the perturbation theory roughly for spherically symmetric, static d-dimensional 
black holes (d > 3), with mass M, charge Q and background cosmological constant A, and the computation 
of QNMs and QN frequencies. We refer the reader to [17] for more detail. 

For a massless, uncharged, scalar field, after a harmonic decomposition of the scalar field as 



- 3 - 



\P = J2i m r 2 V^( r > t) (Oi), where the 8i are the (d — 2) angles and the Ye m (9i) are the d-dimensional 
spherical harmonics, and a Fourier decomposition of the scalar field tpi(r, t) = ^>(r)e lu)t , the wave equation 
can be decoupled Schrodinger-like equation 

where r* is tortoise coordinate defined as dr* = -jfr and V(r*) is the potential, both determined from 
the function /(r) in the background metric. The potential V(r*) depends on the background space-time 
metric and the perturbative type (appendix [A]). For QNMs we need some boundary conditions, so that 

^(r*) ~ e iuJV * as r* -oo, 
^{r*) ~ e~ tujr * as r* -> +00. 

Using this boundary conditions and the monodromy technique, we shall show how to calculate the 
asymptotic QN frequencies and their first-order correction in all static, spherically symmetric black hole 
spacetimes (including asymptotically flat spacetimes and non-asymptotically flat spacetimes). As an ex- 
ample, we may pay more attention to the case of d = 4. For some cases, we shall list some corrected QN 
frequencies, so that we can show it is a reasonable correction by comparing with the numerical results. 

2.1 The Schwarzschild Case 

Although the perturbative calculation for the Schwarzschild solution in 4-dimension and higher dimension 
have been discussed in [22] and [27|, respectively, we first review it roughly for completeness. 



For Schwarzschild black hole, we have 

t( \ 1 2m 

with the roots 

-3 



exp ( -j— -n ) , n = 0, ,d-A. 



(2m) d ~ 

The radial wave equation for gravitational perturbations in the black-hole background can be written as 

d 2 m{u) 



dr 2 



+ V[r(r*)]*(r*) = uj 2 ^(r*) 



in the complex r-plane. 

As mentioned above, the boundary conditions are 

~ e T^ r * as r» ^ ioo, 

assuming Reu> > 0. Then we obtain 

F(r*) ~ 1 as r, -> +00, 

F(r.) ~ e 2iujT * as r* -00, (2.2) 

if we rewrite the QNMs as F(r*) = e tor *^ / (r*). The clockwise monodromy of F(r*) around the r = rn can 
be easily obtained by continuing the coordinate r analytically into the complex plane, i.e., 
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where kn = ^f'{rjj) is the surface gravity at the horizon. 

Near the black hole singularity (r ~ 0), the tortoise coordinate may be expanded as 



dr 



1 r 



d-2 



1 



„2d-5 



+ 



f(r) d-2 2m 2d - 5 (2m) 2 

where f(r) = 1 — and m is the mass of the black hole 1 . When we define z = wr t , the potential 
near the black hole singularity for the three different type perturbations can be expanded , respectively, 
as (appendix |A|) 



V[z] 



"4# 



l-3 2 ~W{j) 



Z \(dr-3)/(d-2) 



UJ 



+ 



(2.3) 



where 



W ST j = 0, 
W(j) = lw sv 3=2, 
Wss j = 0, 

and the explicit expressions of Wst, Wsv and Wss can be found in appendix |A|. Then the Schrodinger-like 
wave equation (2.1) with the potential (|2.3|) wave equation can be depicted as 



Hq + uj d-^HA ^ = 



(2.4) 



where TCq and Tii are defined as 



d 2 

Ho= d^ + 



1-3' 
Az 2 



+ 1 



Hi 



W(j) _a=i 

—Z d-2. 



Taking into account of u — > oo, the zeroth-order wave equation becomes 

H y {0) = 0, 

with general solutions in the form of 

*W=A+J j/2 (z)+A-J_ j/2 (z), 



(2.5) 



where and below J u (x) represents a Bessel function of the first kind. According to the boundary conditions 
(|2.2|) , one can define 

F (0) (z)=F (0) (z) _ e -^/2 F (0) (z)) 

which approaches — e~ ia+ sin 4^ as z — > +oo, where a± = l 1 ^) 71 " . Considering the behavior of the wave 
function as z — ► — oo, we may deduce that 



M(r H ) 



sin 3j7r/2 
sinj7r/2 ' 



1 In fact, one can relate it with the ADM mass M by 



(d - 2) Ad-2 
M — - — — -hz m, 



where An is the area of an unit n-sphere, An 



n + l 
2tv 2 
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leading to the generic expression of (i-dimensional QN frequencies as showed in [17] 

= (2n+l)vri + ln(l + 2cosj7r), 

where Th is the Hawking temperature. 

Next we calculate the first-order correction of the Schwarzschild black hole spacetimes. We first expand 
the wave function to the first order in l/u)( d ~ 3 )/( d-2 ) as 

^ w (d-3)/(d-2) W • 

Then one can rewrite Eq. ( |2.4[ ) as 

W ¥ (1) +Hi^ (0) = 0. (2.6) 

The general solution of Eq. fl2,6|) is 



o 

where C = — . 1 and the wave function are 

sinj7r/2 ' ± 



9>W =»(»)/ <i>!_ 0> Wi<i>i 0) -c*L 0) / ^Mi*"?. (2.7) 



Taking into consideration both the boundary conditions (|2.2j ) and the behavior of wave function F(z) as 
z — * ±oo along the real axis, one may deduce 

, „ , sin3jvr/2 / £_ - K 
M(r H ) . ' 1 + 



sinj7r/2 V w( d - 3 V( d - 2 ) 
leading to the generic expression of (i-dimensional QN frequencies expressed as 

— = (2n + l)?ri + ln(l + 2 cos jvr) + 



T H v 7 v J ; (n + l/2)( rf - 3 V( rf - 2 ) 

where we have defined £± = c-t-te 1 * 14 ^/ 2 — , and 



oc 







And K is defined as 

^ _ e iTr(3+6»y)/4 



sm ( — | — tt) + K- cos( — - — tt) - £ cos( — - — vr) cot — 



where rj = (d - 4) /2(d - 2) and £ = + f_ 



For <i = 4, we obtain the same result as [22 



com = (i - Q 3 ^ + ^ + 1 7 j2 S " (2J?0 * r 2 (i/4)r(i/4 + j/2)r(l/4 - j/2). 

24y2(vr) 3 / 2 sin(3ivr/2) l / M/ J/ M ' J/ 7 

The results above are shown to be dimension dependent and related closely to I and j. It is reasonable 
that these (i-dimensional frequencies would indicate some information about back ground spacetime and 
perturbation types as shown in lowest QNMs. It approaches a constant for the real part of the QN frequency 
once we let n — ► oo as expected in the literature [15, [l6|, |28|. In addition, p^[ showed that the results do 
agree with the result from a WKB analysis p9|], as well as the numerical results |3j for scalar perturbations 
and gravitational perturbations, respective ly. 

2 Throughout this work we have the same definitions for a±, c±, £±, £, and C. 
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2.2 The Reissner Nordstrom Case 

Now we discuss the QNMs of the Reissner-Nordstrom <i-dimensional black hole including first-order cor- 



rections. The calculation for zeroth-order was first done in |L(|, and latter in 17]. In order to state more 
clearly, we start with the zeroth-order calculation. For Reissner-Nordstrom black hole, we have 

„2 



f(r) 



1 



2m 

r d—3 



+ 



,2d- 



with the roots 



in 



m. 



exp 



2vri 
d-3' 



-n 



n 



0,1,- 



,d-4, 



where q is the charge of the black hole 3 . 
Again we define 

F(r.) = e^ r **(r*). 

Then the clockwise monodromy of F(r*) around the outer horizon r = r + can be obtained by continuing 
the coordinate r analytically into the complex plane, i.e., 



M{r + ) 



2ttuj 



where fcj 



i/'(r+) is the surface gravity at the outer horizon. 



Near the black hole singularity (r ~ 0), the tortoise coordinate may be expanded as 



dr 



f(r) 2d -5 



■y j,2d—5 2m 8 

_ + 



3d -8 q 4 



+ 



(2.8) 



One can easily learn from that r* — > oo as q — > 0. In fact, in our procedure for expanding Q2.8j ), we 
have assumed -r <1. As a result, instead of expanding potential to the first order in l/u/ d ~ 3 )' ( 2<J ~ 5 ), we 

must expand it in 1/ [(r^ ) 2d ~ 5 u] ( d ~ 3 )/( M ~ 5 ) _ 

type perturbations can be then expanded, respectively, as (appendix |A|) 

,2 f / „ \ (d-3)/(2d-5) 



After defining z = wr t , the potential for the three different 



V[z) 




J' 



W(j) 



z 



(ro) 



2d- 5 



+ 



where 



W(j) 



W RNT j 
W RNV ± j 

w RNS ± j 



Jt, 

3V±; 
3s±, 



(2.9) 



and the explicit expressions of Wrnt, Wrnv± and W R ^s± can be found in appendix [A]. Then the 
Schrodinger-like wave equation (|2.1|) with the potential (|2.9D can be depicted as 



-\2d-5 



d-3 
' 2d-5 



3 One can relate it with the charge Q by 



(r )- 

2 _ (ri-2)(ri-3)g 2 
V 87rG d 



Wi * = 0, 



(2.10) 
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Figure 1: Stokes line for the Rcissncr Nordstrom black hole, along with the chosen contour for monodromy matching, 
in the case d = 6 ( we refer the reader to fllTf for detail, and a more complete list of figures in dimensions d = 4, 
d = 5, d = 6 and d = 7 ). 



where Ho and Tii are defined as 



y - d " + 



4z 2 



+ 1 



Hi 



4 



Taking into account of u — ► oo, the zeroth-order wave equation becomes 

Wo^ (0) = 0, 

with general solutions in the form of Eq. (|2.5|). According to the boundary conditions 

F (0) (z)=F (0) (z) _ e -^/2 F (0) (z)) 



one can define 



which approaches —e ia+ sin ^ as z — > +oo. This holds at point ^4 in Figure 4 [l]. By going around an arc 

2d- 5 



of angle of Jf^, in complex r-plane (rotating from A to the next branch), z rotates through an angle of 



2ir in z-plane, leading to the wave function 



1 

2 L 



e 2iz f e 3ia+ 



4. i e 5ia + — e ~ i i n /' 2 e 5ia - 



(2.11) 



As one follows the contour around the inner horizon r = r_, the wave function will be of the form 

F<°>(jz - 5) = ^ + Fj 0) (z - 5) - A_f1 0) (z - 5), 



4 A11 figures in this paper are provided by J. Natario and R. Schiappa in Jit| ]. 
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where 

2ujiri uiri 

and = g /'(?"—) is the surface gravity at the inner horizon. Notice that z — 5 < on this branch, one 
can easily obtain 



2 

where we have used 



F^(z -S) = - \e 2l( - z -V (A + e ia + - A_e ia ~) + (A + e~ ia + - A_ e - ia -)] , (2.12) 



Jv{z) = V ^ cos ( z + T + l) ' Z<< " L 



Since on the same branch, we must let F^°\z — 5) = (e 2m z) , and hence it is easily seen from Eqs.(2.11) 
and ( 2.12 ) that we must have 



e 



+ - e -^/2 e 5ia_ = A -ia+ _ A _ e ~i*- _ 



Then we obtain 



a ^ <?it ni J 71 " ,v 3 ~3 sin(3j7r/2) 
A. = 2e 2jd e^~ cos - e J7r ~ 



2 " sin(jvr/2) ' 

A _ = 2e » e *r^ cos £ _ e f siD ( 3 W2) (2 13) 

2 sin(jvr/2) v ; 

Finally, we must rotate to the branch containing point B. This makes z — 5 rotate through an angle of 2tt, 
leading to the wave function 

F^(e 2wi {z - 6)) = 2e 2l5 e~ ia + cos(^) sin(ivr) + e~ ia + sin(^). 



Consequently, we deduce that 



>M = -2e 2i5 [l + cos(jvr)] - (1 + 2cos(jvr)), (2.14) 



which agrees with the result shown in [17|. 

Next we compute the first-order correction of the asymptotic QNMs of the d-dimensional RN black 
hole using the monodromy method. The same reason as mentioned above, we must expand the wave 
function to the first order in 1/ [(rQ~) 2d ~ 5 u;] 3 )/( 2d 5 ) ag 

$ = $(o) + I $(i) 

Then one can rewrite Eq. fl2,10| ) as 

W * (1) +Wi^ (0) = 0. 
The general solution of this equation, as mentioned in the last subsection, is 

= c*s? r *i o) wi*i o) - r ^n^. (2.i5) 



- 9 - 



The behavior as z 3> 1 is found to be 

&± (z) ~ c_± cos(z — a + ) — c + ± cos(z — a_). 

After defining 



* = *< 0) + 



IK 



2d - 5uJ ] (d-3)/(M-5) 



(1) 



-jjvr/2 / j 



M-5, ,1 (d-3)/(2d-5) 



$(°) + 



[(r- 



-\2d-5 



1 (d-3)/(2d-5) 



(1) 



and using boundary condition fl2.2Q , one finds that the wave function -F(z) approaches 



F(z) ~ -e- ia + sin ^ 



IK) 



2</- 



5^1 (d-3)/(2d-5) 



as 2; — » 00. 



The function defined in ( |2.15| ) follows that 



where 7/ = — l/(4d — 10), and G±(z) are even analytic functions of z. By going around an arc of angle of 
2 ^ 5 in complex r-plane (rotating from A to the next branch), z rotates through an angle of 2ir in z-plane, 
leading to the wave function 



F(e 27Ti z) = F ( -°\e 2jTi z) + 



\M-5,.,1 (d-3)/(2d-5) 



-i7r(l+2r;) p(0). 



[(r 

As one follows the contour around the inner horizon r = r_ , the wave function will be of the form 

F(z -5) = A+F+(z -5)- A-F-(z - 5), 

where again 

2iOTii uiri 



f(r_) k- 

and k~ = ^ /'(?■*_) is the surface gravity at the inner horizon. For highly damped QNMs (uj S> 1), 
approximately, we have 

A + = A+, A- = A-. 

Finally, we must rotate to the branch containing point B. This makes z — 5 rotate through an angle of 2tt, 
leading to the wave function 

F(e 2Ti {z - 5)) = e 2t{z - 5) B + i (A + e 3ia + - A.e 3ia ~) 



2r\-Ki 



+ - 



2 [(r -) 



-\2<2-5 



UJ 



I (d-3)/(2d-5) 



[A+C+e 3ia + + A_£_e 3ia - + A^e'^e 



3ia- 
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where B is a constant which is not needed in our calculation because e 2l ^ z ^ — > as z — > — oo. Conse- 
quently, as one uses the explicit expression of A + and in Eq. ( |2.13| ), we deduce that 



F(e 2m (z-5)) = 2e^e- ia cos(V)sin(jvr) 1 



2 ' W ' \ [( r -)2d-5,„l( d - 3 )/( M - 5 ) 



+e sm(— -J 1 



2 \ [(r - )M _ 5a;1 ( d -3)/(2 d -5) 



where 



Finally, we have 



2 -,t .i- 



2e^[l + cos(^)] I 1 + - — == B , I - [l + 2cos(j^)] I 1 + 



Ki = e i7r(??+3/2) [f + sin^vr) - £cos(??7r) cot(jvr) + i£_ cos(r?vr)], 
K 2 = e i7r ( J ?+ 3 / 2 ) [£ + sin(ryvr) - £cos(r?7r) cot(jvr/2) + cos(r?vr)]. 



leading to the generic expression of d-dimensional QN frequencies 

e^+2e ^ [l+cos(j7r)] ( 1 + — ~^\, ll0l1 ^ | +[1+2 cos(jTr)] ( 1 + ^ 5 I = °> 

V [(r -)2^] (d - 3)/(M - 5 V V [(r -)2^] (d - 3)/(M - 5 V 

where and represent the Hawking temperature at the outer and inner horizons (notice that Tjj < 0) . 
For d = 4, jx = js± = 1/3, jy± =5/3, one obtains 



[(r -)3^]^V V [(r ) 



where _ 

Si = l - ± j^[t + ~ i- + V3cot(ivr)], £ 2 = -£- + V3cot(i7r/2)]. 

In the case of g < m, one has — — — * 0, which leads to 

1 H 

W /„ , ,n • , (V3-*)(fi+£2) 

— ^ ~ (2n + l)7Tz + log 5 



T+ v 7 5r„ [(8n + 4)A; + ] 1 /3- 

In the case of q — > m, one has -V « — which leads to 

— r ~ (2n + 1)tti - log 2 + - — — 
T+ V ; & 4r~[(8n + 4)A; + ] 1 /3 

In other cases, approximately, one has 

A + 3e"£ f 1 (V5 - i)g O + 2 ( 1 J = 0, (2.16) 

I [(8n + 4)/c+] 1 /3 I 1 r -[( 8ri + 4 )fc + ]i/3 J ' 1 ' 
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In order to obtain an explicit expression of £ i and £2 , we need the integral 

r(l)r( 



Jo 



M+^+l/3 > 
2 



3/|p ^-/J+5/3 -|p ^+i/+5/3 -j-p ^— i/+5/3 ^ ' 



As a result, we obtain 



;sin(j7r/2) 
!sin(i7r/2) 



J(j/2,j/2)e^/ 2 -J(j/2,-j/2) 
J(-j/2, -j/2)e"^/2 _ J(i/2; _ j/2) 



Immediately, one has 
fl = - 



(V3 - i)^0-) x cos(3 J7 r/2) r(2/3)r(1/6 + j/2)r(1/6 _ i/2)r 2 (1/6)j 



64^2 



cos 2 (j7r/2) 



* 2 = - iV l 2 ^2 ij) • r(2/3)r(i/6 + j/2)r(i/6 - j/2)r 2 (i/6). 



In this way, one can easily obtain the value of QN frequencies in Eq. ( 2.16| ). For tensor type per- 
turbations, Wrnt = 0, so £\ = £2 = 0, leading to the same results as the zeroth-order asymptotic QN 
frequencies. However, it seems unavailable for scalar type perturbations, since in this case we have j — ► 1/3, 
which may induce the integral ^(—1/3, —1/3) approaches infinity. It is interesting to investigate this prob- 
lem in detail. Is there any other methods can avoid this singularity? For vector type perturbations , we 
have j — > 5/3, which lead to 



£1 



0. 



£2 ~ (0.3734 - 0.2156?) (l ± ^9 + M{i + l)q 2 ~ 8q 2 ) (1 - s/l - q 2 )q~ 4/3 , 

with m = 1. From here we find: (1) in order to insure r^[(8n + 4)/c + ] 1 / 3 > 1 as one calculates the QN 
frequencies of the first-order correction, the imaginary part of the frequencies ( or the modes n) needs bigger 
values for a black hole with small charge, since ~ q 2 — > as q — > 0. This confirms the prediction made 
by Neitzke in [3C]: the required n diverges as q — > 0, and the corrections would blow up this divergence; (2) 
just like the case of zeroth-order QN frequencies, the q — * limit of RN corr^ does not yield Schwarzschild 



C07T4, and the same thing happens in the limit of q 



m 



as we shall see in the following; (3) the 



first-order correction to the asymptotic QN frequencies are shown to be dimension dependent and related 
closely to £, j, and the charge q. 



QN frequencies of RN black holes were calculated numerically by Berti and Kokkotas in [31]. They 
found that very highly damped QNMs of RN black holes have an oscillatory behavior as a function of the 
charge. Their results have a good agreement with the zeroth-order analytical formula ( |2.14 ). However, it is 
necessary to perform further checks to our first-order corrected results both in four and higher dimensional 
black hole spacetime. 

2.3 The Extremal Reissner Nordstrom Case 

Now we discuss the QNMs of the extremal Reissner-Nordstrom <i-dimensional black hole including first- 

2 

m 



order corrections. The calculation for zeroth-order was done in |17J, where they found that the q 
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limit of RN QN frequencies does not yield extremal RN QN frequencies. In this case, the outer horizon 
approaches the inner horizon. 
Again we define 



•tf(r*). 



Then the clockwise monodromy of F(r*) around the horizon r = can be obtained by continuing the 
coordinate r analytically into the complex plane, i.e., 



M{r ) = e k o , 



where k, 



. - (d-3) 2 
2{d-2)mM d 



Near the black hole singularity (r ~ 0), the tortoise coordinate may be expanded as 

dr 



f(r) 2d — 5 m 



I r 2d—5 2 r 3d— 8 

+ ~ 5- + 



3d -8 m A 



Again, when we define z = car*, the potential for the three different type perturbations can be expanded, 
respectively, as (appendix [A]) 



z \ (d-3)/(2d-5) 



+ 



(2.17) 



where 



W(j) 



W RNT J=JT, 
W RNV± j=JV±, 
{ W RNS± 3=3 S ±, 



and the explicit expressions of Wf^ NT , W^ NV± and W^ NS ± can be found in appendix [A|. Then the 
Schrodinger-like wave equation (2A) with the potential (|2.17| ) can be depicted as 



H + u 2d=sWi I ^ = 0, 



(2.18) 



where TLq and TL\ are defined as 



d2 

H ° = dz^ + 



1-J 5 

Az 2 



+ 1 



-Hi 



W(j) 3d-7 



Taking into account of u —* oo, the zeroth-order wave equation becomes 

Wo* (0) = 0, 

with general solutions in the form of Eq. fl2.5|) . According to the boundary conditions (|2.2j), one can define 

F(°\z) = F| 0) (z) - e-^F^iz), 

which approaches —e~ ia + sin ^ as z ^ +oo. This holds at point A in Figure^. By going around an arc 
of angle of 2 d-5 m com pl ex r-plane (rotating from A to B), z rotates through an angle of 57r in z-plane, 
leading to the wave function 



F (o) ( 



e 57ri z) 



5jtT 
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contour 




Figure 2: Stokes line for the extremal Reissner-Nordstrom black hole, along with the chosen contour for monodromy 
matching, in the case d — 6 ( we refer the reader to 1 17 for detail, and a more complete list of figures in dimensions 
d = 4, d = 5, d = 6 and d = 7 ). 



Consequently, we deduce that 



2tvu> 8111 



e k 



sin 



jn ' 



leading to the generic expression of (i-dimensional QN frequencies 



2niri + ln(l + 2cosj7r + 2 cos 2jn), 



where T = 2irko. 

Next, we calculate the first-order correction to the asymptotic frequencies. Again we expand the wave 
function to the first order in l/uj( d ~ 3 )/( 2d ~ 5 ) as 

$ = tf(0) _j_ I ^(i) 

Then one can rewrite Eq. ( |2.18| ) as 

H Q ^ {1) +H l ^ ) = 0. 
The general solution of this equation, as mentioned in the last subsection, is 



(2.19) 



The behavior as z 3> 1 is found to be 

&±(z) ~ c-± cos(z — a + ) — c + ± cos(z — a_). 
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After defining 



+ w (rf-3)/(2d-5) * + 



-ijir/2 



1 



*L 0) + 



l 



(i) 



w (d-3)/(2d-5) y ^ - w (d-3)/(2d-5) 

and using boundary condition (|2.2|), one finds that the wave function F(z) approaches 



F(z) e~ ta + sin 



J7T 



1 



w (d-3)/(2d-5) 



as z — ► oo. 



The function VP^ defined in ( |2.19| ) follows that 



^(1) = z l±i/2+r, G±(z); 



± 

where 77 = — l/(4d — 10), and G±(<2) are even analytic functions of z. By going around an arc of angle of 
2%- 5 m complex r-plane, z rotates through an angle of 5tt in z-plane, leading to the wave function 



F(e bm z) = e~' Mz B + i (e 

+ 



2uJ (d-3)/(2d~5) I 



-ijn/2 9io_ 



-5r/7ri — ijir/2 9ia_ 



where B is a constant which is not needed in our calculation because e 2tz — » as — > +00. Finally, we 
have 



e 2 ^ = [1 + 2 cos(jvr) + 2 cos(2jvr)] ( 1 + " A 



3)/(2d-5) / ' 



where 



^ = e (17+10»?)7ri/4 



( . , 3 + 107? 

4+ Sin 7T — 4 cos 



3 + 10?? \ , f 5jn /3 + lOr/ 

^ I COt( — ) + COS I 7T 



leading to the generic expression of d-dimensional QN frequencies 



w corrd 
- = 2nm + In [1 + 2 cosOvr) + 2 cos(2j7r)] + w(d _ 3)/(2d _ 5) 



For d = 4, one obtains 



corr 4 



\/3i - 1 



4^fe 







\/3C+ - v 7 ^- - £cot(5j7r/2) 



(2.20) 



In order to obtain an explicit expression of corr^, we need the integral 

roo p/2MV M+"+V3 

J^fi)^ J dzz- 2 / 3 J u (z)J,(z) = ^ 



r(j)r(S- 



V^4r( ^~^ +5 / 3 )T( At+^+5/3 ^j-i ^— i^+5/3 n 
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As a result, we obtain 



£ + = ^, (J) [j(i/2,i/2)e^/2 _ J(i/2,-i/2) 
8sm(j7r/2) L 



8sin(jvr/2) 
Immediately, one has 



J(-i/2, -i/2)e-^/2 _ J(i/2) _ i/2) 



cos^-cos^-(l-\/3i)cos^ 
8^4^7r 2 X sin sin 52ZE 



co ^4 = - o3/7 f- 9 x — ^ . V r 5j7r — • r(2/3)r(i/6 + i/2)r(i/6 -j/2)r^(i/6). 



In this way, one can easily obtain the explicit expression of corr^ of Eq. (|2.20| ). For tensor and scalar 
type perturbations, W^ NT = W^ NS± = 0, so corr^ = 0, leading to the same results as the zeroth-order 
asymptotic QN frequencies. For vector type perturbations , we have j — ► 5/3. Strangely, this also leads to 
corr<± — ► 0. It is interesting to investigate whether the first-order corrections for any d > 3 extremal RN 
black holes have a same behavior — they all approach zero, independent on the dimension. Another point 
is that in the limit of q 2 — ► m 2 , corr^ of RN black hole does not approaches zero — they have different 
correction terms. 

2.4 The Schwarzschild de Sitter Case 

Now we compute the quasinormal modes of the Schwarzschild dS d-dimensional black hole including first- 
order corrections. We start with the zeroth-order calculation. For Schwarzschild de Sitter(SdS) black hole, 
we have 

with the roots 

r n = r H ,rc,n,rl, ■ ■ • , r d - 4 , r* d _ A , f, 
where A > is the black hole background parameter related to the cosmological constant A by 

A = i(d-l)(d-2)A, 
and r* represents the conjugate of r n . Here we have defined 

(d-4)/2 
t=l 

The clockwise monodromy of 1 f (r* ) around the event horizon r = rn and the cosmological horizon r = rc 
can be obtained by continuing the coordinate r analytically into the complex plane, respectively 

TVUJ 

M(th) = e k n , 

— TVUJ 

M(r c ) = e k c , 
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where ku = hf'( r H) an d kc = hf'(rc) are the surface gravity at the event horizon and the cosmological 
horizon, respectively. 

Near the black hole singularity (r ~ 0), the tortoise coordinate may be expanded as 



f{r) d-2 2m 2d - 5 (2m) 2 



+ 



After defining z = ojr*, the potential for the three different type perturbations can be expanded, respec- 
tively, as (appendix [A|) 



(2.21) 



where 



W(j) = { 



Wsdsr j = jr, 
Wsdsv j=jv, 
Wsdss 3 = js, 



and the explicit expressions of WsdST, WsdSV and Wsdss can be found in appendix [A|. Then the 
Schrodinger-like wave equation (|2.1|) with the potential ( 2.21| ) can be depicted as 



Hq + uj "Wi * = 



where TLq and TL\ are defined as 



ov d2 



I- 3 d 

Az 2 



+ 1 



Hi 



W(j) _a=i 

—Z d-2. 



Obviously, the zeroth-order wave equation can be written as 

H y {0) = 0, 

with general solutions in the form of 

^)(z)=A + ^J ]/2 {z)+A^J_ j/2 {z) 
As one lets z — > +oo, the wave function approaches 
#(o)( 



z ~ 



A +e - ia + + A_e- ia ~) e — + (A+e ia + + A_e ia -) & — 



(2.22) 



This holds at point A in Figure ||. By going around an arc of angle of m complex r-plane (rotating 
from A to B), z rotates through an angle of 3tt in z-plane, leading to the wave function 



¥°\z) ~ (A + e lia+ + A^e lia ') — + (A + e 5ia + + A_e 



5ia~ 
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Im 



Stokes line 




Figure 3: Stokes line for the Schwarzschild de Sitter black hole, along with the chosen contour for monodromy 
matching, in the case of dimension d = 6 (we refer the reader to (ItJ for detail, and a more complete list of figures 
in dimensions d = A, d = 5, d = 6 and d — 7 ). 



Consequently, we have 



A e 7ia+ i A_e 7ia ~ ™ i ™ 

+_ ^e k « + k c = M(r H )M(r c ), (2.23) 

A+e~ ia + +A_e~ ta - 



7VLJ 7TUJ 

c ktr fcr* 



A+e ia + + A_e ia - 



M{t h )M{t c ). (2.24) 



Taking the condition for these equations to have nontrivial solutions (A+, A_) into account, one may easily 
obtain the final results as 

cosh - ™) + (1 + 2a*j«) cosh ( ^ + ^) = 0. (2.25) 
\k c k H J \k H k c J 

As discussed in Ref. [17], the Stokes lines for d = 4 and d = 5 are a bit different from the case discussed in 
the previous calculation for d = 6. On account of the behavior of r ~ oo, they found that the final result 



for d = 4 is the same as Eq. Q2.25 ). However, the result for d = 5 is changed. For r ~ oo, the coefficient 
of e~ lz in the formula of the wave function keeps unchanged, while the coefficient of e lz in there reverses 
sign as one rotates from the branch containing point B to the branch containing point A in the contour. 
Therefore one can complete this calculation by reversing the sign of Eq. ( |2.23j ). In the end, we obtain 

, / 1TUJ iru>\ , v f ITU 1TLd\ 

smh + (1 + 2cosj7r) smh 1 =0. 

\kc kH J \kh kc J 

Next, we calculate the first-order correction to the asymptotic frequencies. Again we expand the wave 
function to the first order in l/^( d_3 )/( (i ~ 2 ) as 

= *(0) + - $(!) 
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Then one can rewrite Eq. ( gj|) as 

H y {1) + Hi^ (0) = 0. 

The general solution of this equation is 

= C< 0) f *L 0) Wi*i 0) - C*i 0) f ^Wi*^. (2.26) 
jo J o 

The behavior as z 3> 1, as mentioned above, is found to be 

*&±(z) ~ c~± cos(z — a + ) — c+± cos(z — a_). 

After defining 

d-3 

a±± = c±±a; d - 2 , u± = 1 ± a_| , 

one finds that the wave function ^f(z) approaches 

~ (A + u + + A-a ) cos(,z — a + ) + (A_u_ — A + a ++ ) cos{z — a_) 

= [{A + u + + A^a^)e~ ia + + (A_u. - A + a ++ ) e - ia ~] ^ + 

[{A + u+ + A_a—)e ia + + (A_u_ - A + a++)e iQ 



2 

as 2 — > oo. 

The function defined in follows that 

where rj = (d — 3)/(c£ — 2), and G±(z) are even analytic functions of z. By going around an arc of angle of 
1X1 complex r-plane, z rotates through an angle of 3tt in z-plane, leading to the wave function 

V(z) ~ [(A + v+ + ,4_a_^e- 3i ^)e 7ia + + {A.v- - A + a ++ e 3ijn )e 7ia -] 6 6 + 

3 -3ij7r\ 5ia + , / 4 _ 4 3yTr\ 5ia_1 c 



+ ^_a__e- 3i ^)e 5ia + + (A„v. - A + a ++ e 3ijn )e 5 



1 2 

as 2 — » —00. Here we have defined v± = e~ 3vm ± We hence have the similar formulae as shown in 

( gjg ) and ( gjg ), i.e., 

^L_L \ —^q;, ■ /\ 1 \ -»a e " c = e ^ ^ ir c ), 2.27 

(A + u + + A-a )e ia + + {A-U- — A + a ++ )e ia - 

(A + v + + J 4„a_e~ 3 ^' 7r )e 5ia + + (A-V- - A + a ++ e 3i ^)e 5ia - ™_«ai _ . , „ x , „ x / 

1 71 Zl ( ia+ T 71 2 >> fa e fcc = e-^M(r H )M(r c ). 2.28 

+ A-a )e ia + + (j4_ii_ — A + a ++ )e ta - 

In this way, one can easily obtain a set of equations with regard to A + and A_ 

s ie 7ia + - s 2 e~^e- ia +^ A + + (s 3 e 7ia - - s 4 e~^ ' e~ ia ~\ A_ = 0, 
■ Sse 5i«+ _ S6e ¥ e »+) A + + f s 7 e 5ia - - s 8 e^e ia -) A_ = 0, 
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where 



IJ7Y q 

si = (v+ - a+ + e~~)e'' nm , s 2 = (u + - a ++ e~), 

s 3 = (v- + a e^)e 3 ' nm , s 4 = («_ + a e -1 ^), 

s 5 = (v + - a ++ e 2 ) e ° ' , s 6 = (u+ - a ++ e 2 ), 
s 7 = (t>_ + a e _i 2 L )e 3,?7r \ s 8 = («_ + a e 1 ^). 



The condition for these equations to have nontrivial solutions (^4+ , A_ ) is then 



s\e 



7iau 



s 5 e 



5i«+ 



s 2 e k c e ia+ 
— fc « e lOL+ 



s 3 e 



7ia- 



S 7 e 



5ia- 



27TU 

s 4 e k c e~ ta - 

2-7TU 

- s 8 e fc H e m - 



0. 



After some algebra, we obtain 



, . TTUJ TTU! . \ . ( TTUJ TTUJ 

Ai cosh h Ai + A% cosh \— + - h A 2 

k c k H ) ~ \k H k c 



(2.29) 



where 



Ai 
A 2 



sis 7 e 2 - S3 s 5 e 2 S4 s 6 e 2 



tJ7T 

s 2 s 8 e~~ 



S3S6e 2~ 



S4S5e _ 2 _ 



s 2 ,s 7 e~~ 



log sis 7 e~ - s 3 s 5 e _ ~ 



jj TV jjlT 

log s 4 s 6 e~ - s 2 s 8 e ~ 



log sis 8 e~ - s 3 s 6 e 2 - log s 4 s 5 e~ - s 2 s 7 e 2 



For zeroth-order asymptotic QN frequencies, one has a±± = 0, and hence u± = 1 and v± = e~ 3vm . As a 
result, formula ( [2.29D reduces to ( 2,25 ). 

As mentioned formerly, the Stokes lines for d = 4 and d = 5 are a bit different from the case discussed 
in the previous calculation for d = 6. Resorting to the behavior of r ~ 00, one finds that the final result 
for d = 4 is shown to be the same as Eq.( 2.2S| ). However, the result for d = 5 is changed. For r ~ 00, the 
coefficient of e~ lz in the formula of the wave function keeps unchanged, while the coefficient of e lz in there 
reverses sign as one rotates from the branch containing point B to the branch containing point A in the 
contour. Therefore one can complete this calculation by reversing the sign of Eq. (ET27D . In the end, we 
obtain 



, , 7TU! TTUJ . \ . . , / ITU TTUJ 

Ai smh h Ai + A 2 smb. — + — + A 2 

kc k H J \k H k c 



0. 



For d = 4 SdS black holes, one has rc 3> rjj as < A <C 1, as a result, one has ku 3> kc- In this 



case, we can neglect kc compared to kn- Taking Eq. (2.25) into account, one finds the asymptotic QN 
frequencies in this case is the solutions of this formula 



e k c + e k c 



0. 
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So, we have ^ = (2n + l)m. However, as one lets A approaches its extremal value, i.e., A — * 1/27 (with 
m = 1), one has ku ~ kc- So, in this limit, we can obtain the asymptotic QN frequencies by solving this 
formula 

e *c + e *c = 0, 

which in turn derives ^ = ^(2ra + where Tc is the Hawking temperature at cosmological horizon. 
This reminds us to conjecture that the frequencies have values between these two extremal values, i.e., 

^~ X (2n+l)™ + Meu;, - < X < 1, (2-30) 
+ C 2 

where x is a parameter closely related to the cosmological constant A. In fact, it would be an interesting 
work to investigate the relationship between these two parameters, by analytical methods or numerical 
ones. Meuj is the real part of the frequency. 

In order to obtain the explicit expressions of A\, A2, Ai and A2, we need the integral 



^ ' 



As a result, we obtain 



si = 1 - 2iE7 (j) ^cos y - sin yj , «2 = 1 - 2«£7 (j) ^cos =y - sin =y 

s 3 = 1 - 2iLT (J) ^cos ^ + sin , s 4 = 1 - 2if7 (j) ^cos ^ + sin 

s 5 = 1-2U (j) ^cos ^ - sin ^\ , s 6 = 1 + 2U {j) ^cos ^ - sin ^ , 

s 7 = 1-2U (j) ^cos ^ + sin ^ , s 8 = 1 + 2U {j) ^cos ^ + sin ^ 

u(j) = W ^ =r 2 (i/4)r(i/4 + j/2)r(i/4 - j/2). 

U; 32vr 3 /2 v /(2n + l)xfcc 
After some algebra, we obtain the explicit expressions of A\, A2, Ai and A2 in terms of o{l/y/n) 

jlT 



where 



A\ ~ 2% sin : 
2 

^,2 — 2^ sin — — 



, , cos377r/2 

1 - 2iU(j) cos jir/2 + 



1 + 2 COS J7T 

A 2 ~ 2E7(j) cosj'tt/2 + 



Ai ~ 0, 

cos 3j7r/2 



1 + 2 COS _77T / 

All type perturbations, including tensor and scalar type perturbations (j — ► + ), and vector type pertur- 



bation (j — > 2) have a same behavior that .4.1 and *4.2 approach zero in the limit. This makes Eq. ( 2.29 ) 
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automatically satisfied. As shown in [17], we first regard j as nonzero, then take the limit as j — > 0. In the 
end, Eq. ( 2.29| ) becomes 



/ 7TUJ iru>\ . i , . , / iru -kuj 8U(j)\ 

cosh U - %) + (3 - 8!f,(j))cosh U + + s) = °- (2 - 31) 

Inserting the expressions of (see Appendix [A]), we can easily obtain the explicit expression of U(j). 

In this way, the first-order correction to asymptotic QN frequencies of SdS black holes can be obtained 
by evaluating Eq. (|2.31 ). The result show that the correction term is closely related to £, A (though the 



parameter x and kc), and of course, n. 



Numerical calculation on this case for very highly damped overtone first appeared in [ 32 1 and later in 



[33]. An analytical formula for four dimensional case was deduced in [R4]. They found that the analytical 
results are in good agreement with the numerical results. However, it is interesting to perform further 
checks to our first-order corrected results both in four and higher dimensional black hole spacetimes. 

2.5 The Reissner Nordstrom de Sitter Case 

Now we compute the quasinormal modes of the RN dS d-dimensional black hole including first-order 
corrections. We start with the zeroth-order calculation. For RN de Sitter black hole, we have 

/( r ) = 1 " ^3 + ^cTe ~ Ar '> 

with the roots 

r n = r + ,r_,r c ,r 1 ,rl,--- ,r d _ 4 ,r* d _ 4 ,r, 
where A(> 0) is the black hole background parameter related to the cosmological constant A by 

A = i(d-l)(d-2)A, 

and r* represents the conjugate of r n . Here we have defined 

/ d-4 



r + +r_+r c + /X r i + 



\ i=i 

The clockwise monodromy of ^(r*) around the outer horizon r = r + and the cosmological horizon r = rc 
can be obtained by continuing the coordinate r analytically into the complex plane, respectively 

7TUJ 

M(r+) = e k + , 
M(r c ) = e k c , 

where k + = ^ f'{ r +) and kc = \ f'(rc) are the surface gravity at the outer horizon and the cosmological 
horizon, respectively. 

Near the black hole singularity (r ~ 0), the tortoise coordinate may be expanded as 



f(r) 2d -5 q 2 3d - 8 q 4 
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In this procedure, we have assumed -r «1, where r 



m — \J m 2 — q 1 ) " " , represents the inner horizon 



2d-5, ,1 (d-3)/(2d-5) 



of the RN black hole. Again we must expand the potential to the first order in 1/ [(r ) 
instead of i/ w (d-3)/(2d-5) _ After defining z = wr t , the potential for the three different type perturbations 
can be expanded, respectively, as (appendix [A|) 



V[z 



, ,2 f / _ \ (d-3)/(2d-5) 



+ 



(2.32) 



where 



W(j) 



W R NdST 3=jr, 
W RNdS v± j = 3v± > 
W RNdS s± j = js± > 



and the explicit expressions of WRNdST, W RNdSV ± and W RNdSS ± can be found in appendix [A|. Then the 
Schrodinger-like wave equation ( |2.1| ) with the potential ( 2.32j ) can be depicted as 



(r ) 2d - 5 - 



d-3 
' 2d-5 



where TCq and TL\ are defined as 



ov d2 



1-f 

4z 2 



+ 1 



-Hi 



Wi ) * = 0, 



W(7') 3d-7 
2d-5. 



(2.33) 



Obviously, the zeroth-order wave equation can be written as 

H y {0) = 0, 

with general solutions in the form of 

¥°Hz)=A + ^J ]/2 (z)+A. 
As one lets z — * +oo, the wave function approaches 



J-j/2(z). 



pt>z>- 

¥°\z) ~ (A + e _ia+ + A_e~ ia -) — + (A+e iQ + +^_e ia ") 



(2.34) 



This holds at point A in Figure |4j. By going around an arc of angle of 2 d-5 * n com pl ex 1 — plane (rotating 
from A to the next branch), z rotates through an angle of 2tt in z-plane, leading to the wave function 

*<°>(z) ~ (A + e 3ia + + A_e 3ia -) ^ + (A + e 5ia + + A_e 5ia ~) 
As one follows the contour around the inner horizon r = r_ , the wave function will be of the form 



_ *) = B + J&=*lj j/2 (z -6)+ B J^lj_, /2 (z - 5), 
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Stokes line 
Im A _ 




Figure 4: Stokes line for the Reissner-Nordstrom de Sitter black hole, along with the chosen contour for monodromy 
matching, in the case of dimension d = 6 ( we refer the reader to JlTj] for detail, and a more complete list of figures 
in dimensions d = A, d = 5, d = 6 and d = 7 ). 



where 

2ujiri lotti 



6 



f'(r-) k. 



and fc_ = ^/'(r_) is the surface gravity at the inner horizon. Notice that z — 5 < on this branch, one 
can easily obtain 

*<®(z - S) ~ (B+e ta+ + B_e ia ~) — — + {B + e~ ia+ + B_e~ ia -) , (2.35) 



where we have used 



/ 2 / w f\ 

u{z) = V^ cos v z + T + V ' *<-!■ 



Since on the same branch, we must let ty(°\z — 5) = ^^(e 2m z), and hence it is easily seen from Eqs.( |2l34D 
and ( 2.35 ) that we have 

A + e 3ia + + A_e 3ia ~ = (B + e ia+ + B.e ia -)e~ l5 , (2.36) 
A + e 5ia + + A_e 5ia ~ = (B +e - ta + + £_ e - iQ - )e i& . (2.37) 

Finally, we should rotate to the branch containing point B. This makes z — 5 rotate through an angle of 
2ir, leading to the wave function 



i{z-S) -i{z-S) 



&°\z - S) ~ {B + e 5ia + + B_e 5ia ~) + (5+e 3ia + + £_e : 
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Consequently, we have 



5io_ 



B+e Ma + +B_e 
A + e~ ia + + A^e- ia - 

B + e 3ia + + B_e 3ia - 



_ e ~iS e Z+ + Zc = M(r+)M(r c ), 



. e iS e kH kc = M(r+)M(r G ). 



(2.38) 



(2.39) 



A + e ia + + A^e 1 

Taking the condition for these equations ( j2.36 -^ 2.39 ) to have nontrivial solutions (A + , A_, B + , U_) into 
account, one can easily obtain the final results 



ttuj ituj\ 



ircu ircu 



A: 4 



k 



c 



( ircu 2iroj 



7TU) 



cosh — +(1 + 2cos(i?r)) cosh — + — +2(l+cos(jvr)) cosh — + - — + — =0. (2.40) 



kcj 



k + k c J 

As discussed in ref jlT]] , the Stokes lines for d = 4 and d = 5 are a bit different from the case discussed 
in the previous calculation for d = 6. On account of the behavior of r ~ do, they found that the final result 
for d = 4 is the same as Eq.( [2.40 ), However, the result for d = 5 is changed. For r ~ oo, the coefficient 
of e~ lz in the formula of the wave function keeps unchanged, while the coefficient of e lz in there reverses 
sign as one rotates from the branch containing point B to the branch containing point A in the contour. 
Therefore one can complete this calculation by reversing the sign of Eq. ( 2.38| ). In the end, we obtain 



1TCU 7UaA 



k c k + J 



7TLO 7TU) 



sinh — + (1 + 2 cos(2vr/5)) sinh — + — + 2(1 + cos(2vr/5)) sinh h - — + — 



1TCU 



2ircu ircu \ 



A-4 



kc J 



0. 



Next, we calculate the first-order correction to the asymptotic frequencies. Again we expand the wave 

2d-5, ,1 (d-3)/(2d-5) 



function to the first order in 1/ [(r Q "cu\ as 

1 



\J/ = \J/(°) _|_ 

Then one can rewrite Eq. fl2,33| ) as 



2d _ 5 , <d-S)/<M-6) 



[('0 

H y {1) +n 1 ^ {0) = o. 



The general solution of this equation is 

Jo Jo 
The behavior as z S> 1, as mentioned above, is found to be 

ty±\z) ~ C-± COs(z — Q + ) — C + ± COs(z — Q 

After defining 

\2d-5 



(0) 
± • 



(2.41) 
(2.42) 



a±± = c±± 



(To) 4 



10 



d-3 
2d-5 



u± 



l±a+_, 



one finds that the wave function ^f(z) approaches 

*{z) = ^°\z) + [(r ) 2d - 5 u;Y M *«(*) 

~ (A + u + + A-a ) cos(z — a + ) + (A_u_ + A + a ++ ) cos{z — q_) 

= [(A + u + + A^a^)e- ia + + (A_«_ - A + a ++ ) e - ia ~] ^ + 
[(A + u + + A_a^)e ia + + (A_u_ - A + a++)e ia -] ^ 
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as z — > oo. 

The function ty± defined in ( |2.42| ) follows that 

= z 1± ^2+ ??G±(z)) 

where rj = —l/(4d — 10), and G±(z) are even analytic functions of z. By going around an arc of angle of 
2 J^ 5 in complex r-plane (rotating from A to the next branch), z rotates through an angle of 2ir in z-plane, 
leading to the wave function 

~ - ,4_a__e- 2 ^)e 3ia + + {A_v + + ^+a ++ e 2 ^ 7r )e 3ia -] 6 6 + 



2 

- J 4_a__e- 2i ^)e 5iQ + + + A + a ++ e 2ijn )e 5ia ~] "-—^ — . (2.43) 



as z — > — oo. Here we have defined v± = e~ 2r)m db a_| 

As one follows the contour around the inner horizon r = r_ , the wave function will be of the form 

_ d-3 

q?( z -5) = ¥°\z -5)+ Ur ) 2d - 5 J 2d_s *W(z - 5), 

where 

2am Ldiri 

and = |/'(r_) is the surface gravity at the inner horizon. Notice that z — <5 < on this branch. As one 
lets z — 5 — > — oo, the wave function becomes 

¥(* - <5) ~ [{B + u + + £_a__)e iQ + + (5_u_ - B+a ++ )e fa -] ^— ^ + 

[{B + u + + J B_a__)e- ia + + (£_u_ - B + a ++ )e~ ia -] 6 ^ - , (2.44) 

where we used 

= V^ cos v z + T + iJ ' Z «- L 



Again, since on the same branch, we must let ^f(z — 5) = \&(e z). Then it is easily seen from Eqs.( 2.43 ) 
and (2.44) that we must have 

[(A + v- - A_a—e- 2ij7T )e 3ia + + (A_v + + A + a ++ e 2ij7r )e 3ia -] e 2r > ni = 

[(B + u + + B_a__)e ia + + (B_u_ - B + a ++ )e ia -] e~ iS \ (2.45) 

[(A+v- - A^a—e- 2ijw )e 5ia + + (A_v + + A + a ++ e 2ij7r )e 5ia -] e 2r > ni = 

[(B + u + + J B_a__)e- iQ + + (£_u_ - B + a ++ )e~ ia -] e tS . (2.46) 

Finally, we should rotate to the branch containing point B. This make z — 5 rotate through an angle of 
2ir, leading to the wave function 

p 2riTri p i(z-6) 

$(z - 8) ~ [(B+v- - B^a^e- 2l ^)e 5ia + + (B_v+ + B + a ++ e 2i ^)e 5ia -] - 



[{B + v_ - B_a__e- 2ijw )e 3ia + + (B_v + + B + a ++ e 2ijn )e 3ia ' 



2 

e 27]iri e -i(z-S) 
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Consequently, we have 
{B+v. 



B_a__e- 2i ^)e 5ia + + (B_v + + B + a ++ e 2i ^)e 5ia - ^^^^ 



(A + u + + A_a )e~ ia + + (A_u_ - A + a ++ )e~ ic 



M(r + )M(r c ),(2A7) 



(B+v- - B-a- 



-2ij7r\ 3ia+ 



+ (B..V+ + B + a ++ e^)e 3ia - w « -g-g 



e e k h 



fee = A4(r+)A4(rc-).(2.48) 



(yl + n + + A_a )e iQ + + (A-U- - A + a ++ )e ia - 

The condition for these equations ( |2.45| ^ 2.48 ) to have nontrivial solutions (A+, A_, B + , B_) is then 



s 2 e 



s 3 e 

ia 



3ia+ 
5ice+ 



2-KU J 



e K c s 4 e 



s 7 e 

ia 



3ia- 




-a 8 e i(a - 


--5) 


5ia- 


_ S2e -i(a+-<5) 


-a4e _i ( a - 


---5) 


-g k C 


_ S3 e*(5«+-<5) 




-5) 




_ Sie i(3a++6) 


_ S7 g*(3a- 


+<5) 



0. 



where 



si 


= («_ 


+ a + 


Z17T o 

+ e ^r)e 2 ^\ 


S 2 = 




SJ7T 

e 2 , 


S3 


= («_ 


+ a + 




S4 = M. 


_ +a_ 


_e 2 


S5 


= («+ 


— a_ 




S 6 = U + 


- «++ 


e 2 , 


•57 


= («+ 


— a_ 




Sg = U. 


_ +a_ 


_e 2 . 



After some algebra, we obtain 



A\ cosh 
where 

Ax 

A 2 
A 3 



7TLO TTLU 

k+ k c 



. , , / 7TCJ 27R<J 7TO; . \ , , ,7R<J 7TCJ . 

Ai ) + A 2 cosh ( — + — + — + A 2 J + A 3 cosh( — + — ) = 0, 



(2.49) 



[{2 Sl s 3 s 5 s 7 - s\s^ - s\s\e-^) (2s 2Si s 6 s 8 - s 2 s 2 e^ - si^e"^)] * , 

[(2s 3 s 5 s 6 s 8 - s 2 s 2 e 2 ^ - 44e~^) {2s 1 s a a 4t s 7 - s 2 s 2 e 2 ^ - s 2 s 2 e' 2i ^)] * , 
S3S 4 s 6 s 7 e l: " r + sis 2 s 5 s 8 e~ lj7T - sis 3 S4S 8 e 2lj7T - s 2 s 5 s 6 s 7 e~ 2ljn , 



Ai = - [log (2 SlS3 s 5 s 7 - s 2 s 2 e^ - slsle-^) - log {2s 2 s iS& s 8 - s 2 s 2 e 1 ^ - s 2 s 2 e -^)] , 



A 2 



[log (2 S3 s 5S6 S8 - sjsje 2 ^ - s 2 s 2 e~ 2i ^) - log {2 Sl s 2 s 4 s 7 - s\s\e 2 ^ - s^e" 2 ^)] . 



For zeroth-order asymptotic QN frequencies, one has a±± = 0, and hence u± = 1 and v± = e 2rim . As a 
result, formula ( |2.49|) reduces to ( 2.40 ). 



As discussed in Ref . |T7|] , the Stokes lines for d = 4 and d = 5 are a bit different from the case discussed 
in the previous calculation for d = 6. On account of the behavior of r ~ do, they found that the final result 
for d = 4 is the same as Eq. Q2.49 ). However, the result for d = 5 is changed. For r ~ oo, the coefficient 
of e~ lz in the formula of the wave function keeps unchanged, while the coefficient of e lz in there reverses 
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sign as one rotates from the branch containing point B to the branch containing point A in the contour. 
Therefore, one can complete this calculation by reversing the sign of Eq. (2.47). In the end, we obtain 



, / TTUJ TTUJ . \ , ( TTUJ 2TTUJ TTUJ . \ , , TTUJ TTUJ . 

A ' smh {To ~ T + + Al ) + M smh ( ~ + + IT + k~ c + A2 ) + ^ 3 smh( ~ + + ~ c ) = °> 

For d = 4 RNdS black holes, one has rc r+ > r_ as < A <C 1. As a result, one has — > fc + 3> &c. 
In this case, we can neglect l/fc+ and —1/k- compared to l/fcc- Taking Eq. fl2.40| ) into account, one finds 
the asymptotic QN frequencies in this limit is the solutions of this formula 

TTUJ TTUJ 

e k c _|_ g k c — Q 

So, we have^- = (2n + l)vri, where Tq is the Hawking temperature at cosmological horizon. However, as 
one lets A approaches its extremal value, i.e., 

__ 2 + 2^9 - 8g 2 
~* (3 + 79^8^)3' 

with m = 1. |3l| showed that in this extremal case we have kn ~ fee — » 0, while does not approach zero. 
So, in this limit, one can neglect the term with regard to 1/k-, and obtain the asymptotic QN frequencies 
by solving this formula 

e <=c + e ~ *c = 0, 

which in turn derives ^ = ^(2n + l)7ri. This reminds us to conjecture that the frequencies have values 
between these two extremal values, i.e., 

^-~ X (2n + l)™ + Meu;, \<X<^ (2-50) 

where x is a parameter closely related to the cosmological constant A. In fact, it would be an interesting 
work to investigate the relationship of these two parameter, by analytical methods or numerical ones. Me u 
is the real part of the frequency. 

In order to obtain the explicit expressions of A\, A2, Aj and A2, we need the integral 

dzz-^ Jjz) Jjz) = ^ [ 2 > 



As a result, we obtain 

ai = 1 - 2U(j) sin ^±^1, S2 = 1 + (i _ v^ )C /(,-) sin £ + 

S3 = l + (l + V3i)C/(j) sin (5 + 3j> , S4 = 1 + (1 - V^i)U(j) sin (5 ~ 3j> , 

fa o 

S5 = 1 + (1 + VZi)U(j) sin (Szi^L, S6 = 1 + 2U(j) sin ^±M^, 

fa O 

sr = 1 - 2C/(i) sin S8 = 1 + 21TC-) sin 

b fa 



r , , , 
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where 



U(j) 



^(j)r(2/3)r 2 (i/6)r(i/6 + j/2)r(i/6 - j/2) 



16ir 2 rQy/(4n + 2) X kc 

After some algebra, we obtain the explicit expressions of Ai, A2, Ai and A2 in terms of o(l/y/n) 
2(1 - cos jvr)[l + (1 + V3i)U(j) cos jir/2], 



A 2 
A 3 

A 2 



2(1 - cos 2jvr) <^ 1 + 



cosjvr/2 - (3/2 + y/&) 



COS JIT 

cos jvr/2 



U(j) 



2(cos jvr — cos 2jvr) 1 1 + 
(3\/3i - 5)C/(j) cosjvr/2, 



4 cosjvr/2 + 2\/3', 



2 cos 2 jvr — cos jvr 
1 + 2 cos jvr 



U(j) 



U(j) 



(2 + \/3i)cosjvr/2 + 



y/3i cos jvr 
2 cosjvr/2 



For vector type perturbations, we have j = 5/3 = 2 — 1/3. Therefore, as one inserts j = 5/3 into 
Ai, A2, A3, Ai and A2, we have exactly the same results as in the j = 1/3 case(tensor or scalar type 
perturbations), and consequently Eq. fl2.49Q becomes 



cosh (j^- - j£ + (5^3/2 - 3i)U(j)j + 

+ [2 + (5v / 3 + 3i)C/(j)]cosh(— + — ) = 

k + k c 



(3y/3 + 15i)U(j) 



cosh (_ + -^— + _ + ( v / 3 + 2i)U(j) 



(2.51) 



for three type perturbations. Inserting the expressions of W(j) (see Appendix [A]), we can easily obtain the 
explicit expression of U(j). In this way, the first-order correction to asymptotic QN frequencies of RNdS 
black holes can be obtained by evaluating Eq. fl2.5ip . The results show that the correction term is closely 
related to t, q, A (though the parameter x and kc), and of course, n. 

Numerical calculation of this case has not been performed in the previous literature. Consequently, 
it is interesting to perform numerical checks to our first-order corrected results both in four and higher 
dimensional black hole spacetime. 

2.6 The Schwarzschild Anti de Sitter Case 



This case was studied analytically in |23j. Here we list their results for completeness. Considering the 
behavior of the black hole singularity (r ~ 0) and large r, they obtained the QN frequencies for all type of 
perturbations, 



VT y . \ 1 % 

-(2 + J + 7oo) - tan" 1 rr- 

4 y J 1 + 2e % ^ 



+ nvr, 



where = d — 1, d — 3 and d — 5 for tensor, vector and scalar perturbations, respectively, and z = ur* is 
the integration constant of the tortoise coordinate (refer the reader to appendix [A]). 
By expanding the wave function to the first order in l/^/ 1 ^ 3 )/^ -2 ) as 



\fr = vEr(°) + 



1 



w (d-3)/(d-2) 
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Ref. [23] showed that the first-order correction to the QN frequency is 



z = ur* = j(2 + j + joo) + + rnr - ^ 6ib - 2ie ln d ~ 2 b - 9ai + a 2 + e lnd ~ 2 (ai - a 2 ) 

where 

ai = ai(oo), 02 = 02(00), 6 = fe(oo), 
and ai(z), 02 (2), are defined as 

ai(z) = — ; v ' 



u ^-2 / dz'z' d-2j j/2 (z')J j/2 (z'), 



a 2 (z) = I dz'z'-^N l/2 (z')N l/2 (z'), 

b(z) = ^^oT^i / dz'z'~^J j/2 (z')N j/2 (z'). 



8 



W(j) in these formulae is defined as 



W(j) 



'WsAdST 3=3T, 
WsAdSV j = jv, 
WsAdSS j=JS, 



and the explicit expressions of WsAdST, WsAdSV and WsAdSS can be found in appendix [A|. 
For d = 4, one obtains 



7T 

= ^ (2 + j + joo) + + corr 4 , 



(2.52) 



where 



corr4 



In 2 _ (l + i)W(j) fK 
2i 128tt 2 ' V n 

and we have used the integral 



7T 1 

oot-(--j)-3 



dzz- l l 2 J u {z)J^{z) 



( cos ^ _ sin jV2)r 2 (l/4)r(l/4 + j/2)r(l/4 - j/2), 



r(i)r( ^ 2 +1/2 ; 



In this way, one can easily obtain the explicit expression of corr^ of Eq. (|2.52[ ) . ]^3| have given the explicit 
expressions for tensor, vector and scalar type perturbations in d = 4 case. Moreover, they showed that 
their results are in good agreement with the numerical results performed in [36]. 

2.7 The Reissner Nordstrom Anti de Sitter Case 

We now compute analytically the quasinormal modes of the RN AdS d-dimensional black hole including 
first-order corrections. Previous work on this case can be found in 37, 39]. Here we start with the 
zeroth-order calculation. For Schwarzschild Anti-de Sitter black hole, we have 



if \ 1 2m X 2 



(2.53) 
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with the roots 

r n = r +) r_ri,r*, • • • ,r d _ 3l r* d _ 3 , 
where A(< 0) is the black hole background parameter related to the cosmological constant A by 

A = i(d-l)(d-2)A, 

and r* represents the conjugate of r n . Near the black hole singularity (r ~ 0), the tortoise coordinate may 
be expanded as 

dr 1 r 2d ~ 5 2m r 3d ~ 8 



+ 



f(r) 2d -5 q 2 3d - 8 q 4 
In this procedure, we have assumed where 



+ 



m — \J m 2 — q 2 



i 

d-3 



represents the inner horizon 

2d-5, ,1 (d-3)/(2d-5) 



of the RN black hole. Again we must expand the potential to the first order in 1/ [(rQ~) 2d ~ 5 w 
instead of l/cj( d ~ 3 )/( M ~ 5 ) . After defining z = wr t , the potential for the three different type perturbations 
can be expanded, respectively, as (appendix [A|) 



V[z]~-^il-f-W(j) 



(d-3)/(d-2) 



+ 



(2.54) 



where 



WRNAdST j = JT, 
W(j) = I W RNAdsv ± j = jy±, 

> W RNA dss± j = 3s± > 

and the explicit expressions of WrnacIST, ^RNAdSV ± an d ^RNAdSS ± can be found in appendix [A|. Then 
the Schrodinger-like wave equation (2.1) with the potential fl2,54|) can be depicted as 



n + 



d-3 

d-2 



Hi ] * = 0, 



where TLq and TL\ are defined as 



d 2 



Ho = dT 2 + 



Az 2 



+ 1 



Hi 



W(j) 3d-7 



Obviously, the zeroth-order wave equation can be written as 

n ^ {0) = 0, 

with solutions in the form of Bessel functions 



As one lets z — > +oo, it behaves as 

¥<°)(*) ~ ( 



^j j/2 (z)+AJ^N j/2 (z) 



~ ia + - iAe~ ia +) — + (e ia + + iAe ta +) - 



(2.55) 
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Figure 5: Stokes line for the Reissner-Nordstrom Anti-de Sitter black hole, along with the chosen contour for 
monodromy matching, in the case of dimension d = 6 ( we refer the reader to ]17| for detail, and a more complete 
list of figures in dimensions d = 4, d = 5, d = 6 and d = 7 ). 



This holds at point £> in Figure ||. Appendix [A| shows that the potential for large r behaves as 

v = joo ~ 1 +■■■ 

4(z - z) 2 

where joo = d— l,d — 3 and d — 5 for tensor, vector and scalar perturbations, respectively, and z is the 
integration constant of the tortoise coordinate 

f°° dr 

~ z = u,f ' =w L WY 

Then the Schrodinger-like wave equation (|2.1|) with this potential can be depicted as 



d 2 *oo(z) 
dz 2 



+ 



1 - i 2 

4(z - zf 



+ 1 



Vooiz) = 0. 



On account of the boundary condition ^^(z) — ► as r — ► oo, one can obtain the asymptotic behavior of 
the wave function 

T / n s r 



e ^ e -i(*+/9) + e -iz e i{z-l3) 



(2.56) 



where /3 = ^(1 + joo)- Since on the same branch, we must let ^°\z) = ^^(z), then one finds 

A = tan(z — (3 — a+). 



(2.57) 



As one rotates from the branch containing point B to the branch containing point A in the contour, through 



an angle of 



2d- 5 



Z ~ 



(2d-b)q 



7 rotates through an angle of — 7r, leading to the wave function 



p-ia+ p iz 

[I - iA) - + [1 - »A (1 + 2e^)] 
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as z — > — oo. 



From the boundary conditions (2.2), we have 

1 - %A (1 + 2e^ ff ) = 0. 
Consequently, taking Eq.( |2.57j ) into account, one obtains the QN frequencies 

7T , _i i 
2 = — (2 + 7 + 700) — tan - — h n7r 

Next, we calculate the first-order correction to the asymptotic frequencies. Again we expand the wave 



function to the first order in 1/ [(r 



2d-5, ,1 (d-3)/(2d-5) 



as 
1 



r 2d _ 5 , ( d -3)/(2 d -5) 



Then one can rewrite Eq. ( |2.55| ) as 
The general solution of this equation is 

TT^/zJj/ 2 (z 



[( 



H ^ {1) +Hi* (0) = 0. 



7T^N j/2 (z) 



dz'V?,L /2 (z')H 1 ^ 0) (z'). 



The behavior as z — > +00 is found to be 



Here we defined 
where 



^ iJ (z) ~ (A ie - ia+ - iA 2 e~ ia+ ) & — + {A ie ia+ + a 2 e ia +) ^ 



Ai = l-b-Ad 2 , A 2 = (1 + 6)A + ai, 



01 = ai(oo), 02 = 02(00), 6 = 6(00), 
and ai(^), 02(2), K z ) are defined as 



<l2{z) 



b(z) 



irW(j) 



irW(j) 



7TW(j) 



) 2d -\ 

(r -) 2d - 5 ^ 
(r -) 2d - 5 ^' 



d-3 r z 
2d-5 



dz z 2d - 5 Jj/ 2 (z)J j/2 (z), 

dz'z'-^N j/2 (z')N j/2 (z'), 

dz'z'-^J j/2 {z')N j/2 {z'). 



This holds at point B in Figure [5|. For the same reason have been discussed above for the zeroth-order 
case, we should let ^f(z) = ^^(z). Combining Eq. (|2.56|), one finds that 



A 2 = A\ t&n(z — (3 — a+) 



as z — > 00. 
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Again, as one rotates from the branch containing point B to the branch containing point A in the 



contour, through an angle of 2^=51 
function (as z — * —00) 



z ~ 



(2cZ-5)g 2 



rotates through an angle of —it, leading to the wave 



$(0)(V) „ + |[i _ iA (1 + 2e^)] e- 3iQ + + e - ia +e- in ^s ( Bl - iB 2 ) 



2 ' 



where 



Bi 



Bo 



-2 cos ^ [e- y,r oi + (A - e" 2ia +) [b- i (1 + e-tf*) o x ] 



e^oa - 4 cos 2 — a x - 2i (l + e ijn ) b 



-2ia- 



[ e -«^oi + A [6 - i (1 + e~ ijv ) oj] , 



and C is a constant that can be easily calculated but is not needed here. 
Taking the boundary conditions fl2.2j ) into account, we have 

[1 - %A (1 + 2e ijn )] e~ 3ia + + e^+e'^^s (Bi - i5 2 ) = 0. 
Consequently, combining Eq. fl2.57l) , one obtains the QN frequencies 



z = ur* 
For d = 4, one obtains 

where 



n , . In 2 1 r 

I (2+J+^) + — +nvr-- 



6ib — 2ie in 2d ~ 5 b — 9ai + a 2 + e lw 2d ~ 5 (ai — a 2 ) 



7T 

w n f» = -(2 + j + j c 



tan 



1 + 2eW* 



+ n7r + corr4, 



(2.58) 



corr 4 = — ln(cosj7r/2) — _ — W 

2i 10247r J r V n7r l 

•r(2/3)r 2 (i/6)r(i/6 + j/2)r(i/6 - j/2) 

and we have used the integral 



(5\/3 + l)sinj7r/ 2 - (5 - V3) cosjvr/2 



/"CO 

/ dzz'^J^J^z) 
■> 



r(|)r( 



M+^+1/3 n 
2 / 



For d = 4, the roots of the formula /(r) = can be depicted by r + (the black hole outer horizon), r_ (the 
black hole inner horizon), and two complex roots, n and its complex conjugate r*. As discussed in p^| , it 
is convenient to set A = —1. From ( |2.53| ), the roots are 



ri 



(5 + y/-2-g 2 + 4m/g), 



-(5- V / -2-5 2 + 4m/ 5 ), 



7t(-5+ \/-2- <? 2 -4m/g), rj = -(-5 - ^-2 - 5 2 - 4m/ 5 ), 
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with g = J— | + ^7=-, where 



ft± = ( 2 - 72g 2 + 108m 2 ± 4(1 + 12g 2 ) 3 + (2 - 72g 2 + 108m 2 ) : 



The integration constant in the tortoise coordinate can be then solved as 
dr 



Sir) 

^,2 

o In 5 In -= — In — 

3r_ + r_ — r + rirj r + 3rf + r\ — r + r_rj r + 3^ d + r* — r+r_n r + 



Obviously, f* is the formula (105) in [23| as g = (or the Schwarzschild limit). 

In this way, one can easily obtain the value of QN frequencies in Eq. fl2.58|) . For tensor type pertur- 
bations, WrnAcLST = 0, so corr^ = — i-^-. However, it seems unavailable for scalar type perturbations, 
since in this case we have j — > 1/3, which may induce the integral dzz~ 2 / 3 J_ 1 / 3 (z) J_i/ 3 (z) approaches 
infinity. It is interesting to investigate this problem. Is there any other methods can avoid this singularity? 
For vector type perturbations , we have j — > 5/3. This leads to 



corn = + -^(1 + ^)^^ 3/2^.^ 

4 1024tito V revr w h 

where 

W RNAdSV ± = \{m± v / 9m 2 +4£(£ + l)g 2 -8g 2 ) (1 - sfl - q 2 )(3q 2 )- 2 / 3 

From this formula we find: (1) in order to insure rQ~(4n + 2) 1 / 3 > 1 as one calculates the QN frequencies 
of the first-order correction, the imaginary part of the frequencies ( or the modes re) needs bigger values 
for a black hole with small charge, since r^ ~ q 2 — > as q — > 0. This confirms the prediction made by 



Neitzke in [3C]: the required n diverges as q — > 0, and the corrections would blow up this divergence; 
(2) the first-order correction to the asymptotic QN frequencies are shown to be dimension dependent and 
related closely to m, £, j, and the charge q. 



QN frequencies of RN Ads black holes were first calculated numerically by Berti and Kokkotas in [40], 



in the case of d = 4. [41] latter performed an extensive numerical study of QN frequencies for massless 
scalar fields in d = 4 RN AdS spacetime. They find for the higher modes, both Mew and Imw increase with 
n. As q increases, Imw increases faster than Re a; does, with n. [17] matched these predictions to their 



zeroth-order analytical formula of the asymptotic QN frequency. Their results are complete agreement with 
the numerical calculations of [41]. In our first-order corrected formula, the frequencies can be obtained by 
subtracting from their results. This leads, obviously, to the same gap of uj as [17] does. Although 



more complete numerical results are needed to check our analytical results in higher dimensional cases, our 
results are in good agreement with the numerical results in the particular case of d = 4. 

3. Conclusions and Future Directions 

In this paper we studied analytically quasinormal modes in a wide variety of black hole spacetimes, including 
d-dimensional asymptotically flat spacetimes and non-asymptotically flat spacetimes. We extend the 
procedure of flTJ to include first-order corrections to analytical expressions for QN frequencies by making 
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use of a monodromy technique, which was first introduced in [16| for zeroth-order approximation, and later 
extended to first-order for Schwarzschild black hole in [22]. The calculation performed in this paper show 
that systematic expansions for uncharged black holes include different corrections with the one for charged 
black holes. In other words, d-dimensional uncharged black holes have an expansion including corrections 
in l/u)( d ~ 2 '' ( d ~ 3 ), while charged ones have an expansion in l/u/ d ~ 3 )/( M-5 ) . This difference makes them 
have a different n-dependence relation in the first-order correction formulae. 

The method applied above in calculating the first-order corrections of QN frequencies seems to be 
unavailable for black holes with small charge, since the required n diverges as the charge q — > 0, and worse, 
the corrections would blow this divergence up. More extensive investigation is needed for this problem. A 
very recent work on RN black holes appeared in [42], where they discussed a possible way to avoid this 
divergence for zeroth-order case in the small charge limit. 

Some remarks are due for the particular case of d = 4 dimensional spacetimes, as in these cases we 
discussed more deeply. For Schwarzschild case, we obtained a j and ^-dependence correction term, which 
was proved by numerical results. For RN case, a puzzle appears when we apply our method to calculate 
the first-order corrections of the scalar type perturbations, since the correction in this case approaches 
infinity. In fact, this problem also exists in other four dimensional charged black holes. Strangely enough, 
extremal RN black hole can avoid this problem. Moreover, it seems that the correction term equals to zero, 
independently of the type perturbations. An investigation on whether the first-order corrections for any 
d > 3 extremal RN black holes have this behavior, should be further performed. Another point worthy of 
being mentioned about RN case is that the extremal limits (q — > 0, or q — > m) of the RN QN frequencies 
do not yield the right QN frequencies in their extremal cases (Schwarzschild QN frequencies, or extremal 
RN QN frequencies). Our results show that the same thing happens in the first-order corrections. |43| 
found that the limit is a singular one since it involves topology change at the level of the contours in the 
complex plane. 

On the applications to quantum gravity [17] indicate that the In 3 in d = 4 Schwarzschild seems to 
be nothing but some numerical coincidences. Our results support this argument. However, a recent paper 
on quantization of charged black holes [44] found that other than the Schwarzschild case, the asymptotic 
resonance corresponds to a fundamental area unit AA = 4Mn2. According to Dreyer's conjecture [13], 
one can immediately set the spin network unit as 1/2, i.e., jmin = 1/2- This differs from j m in = 1 for 
Schwarzschild case. From the LQG point of view, this result support the claims [14] that the gauge group 
of LQG should be SU(2) in spite of the In 3 for Schwarzschild. Is there any other explanation of that? Or 
more generally, can we find other similar proof to support Dreyer's conjecture for other spacetimes? 

There are some other possible directions for further study: (1) perform some numerical calculations 
to check our analytical solutions done above; (2) make sure if it is possible to have such an easy relation 
between the asymptotic QN frequencies and the overtone modes n as suggested in the asymptotically dS 
cases d2~30l) and ( gjSOD , or more deeply, find the possible relation between x and A and a possible explanation 
of that; (3) find an explanation of our correction results on the physical side. 
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A. The Expanded Master Equation Potentials 



In this appendix we shall deduce the expressions of the potentials in Schrodinger-like equation (p7T 



Natario and Schiappa [17] have deduced the value of the potentials in the region of large r. For complete- 
ness, we may list them in this appendix. Our work focuses on the region of the origin, which plays an 
key role in our procedure to calculate the first-order correction of the asymptotic QN frequencies using 
monodromy method. 

We start with a set of potentials determined by the type of perturbations and can be deduced from 



a set of master equations which were derived by Ishibashi and Kodama in [|24| , 25, |26| (can be denoted as 
the Ishibashi-Kodama master equations). Making use of these potentials, we make an exhaustive study of 
the first-order values of all these potentials in terms of the tortoise coordinate. The zeroth-order values of 
these potentials can be found in |17]| . For tensor, vector and scalar perturbations, we obtain, respectively, 

vrw - m (i(w> + c - 2 ^ 3 m + ZzMm . (a,, 



V M ft J W + d-3) (d-2)(d-4)/(r) (d-l)m (d - 2) 2 q 2 A \ 



where 



A = ^(d - l) 2 (d - 3) 2 m 2 + 2 (d - 2) (d - 3) ^ (i + d - 3) - (d - 2) )g 2 . 



(A.3) 



The expressions of C/± and H± can be found in [17|. In these expressions, the ^+ equation represents the 



electromagnetic mode and the equation represents the gravitational mode. 

The Schwarzschild Case: In this case, the black hole is uncharged, so one has Vy+ = Vy- and Vg+ = Vg- . 
For convenience, we denote them by Vy and V$, respectively. The background space-time metric is 

tt \ 1 2m 

Near the black hole singularity (r ~ 0), the tortoise coordinate can be expanded as 



f(r) d- 2 2m 2d - 5 (2m) 2 



+ 



where we have kept the second term in the expansion of r and have chosen the integration constant so 
that r* = at r = 0. Using this formula, one may obtain an approximate expression of r in terms of the 
tortoise coordinate r* 

r ~ [-2m(d - 2)r«]^ + 2) [* • (A.5) 
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As one applies ( |A.5| ) to the three different potentials (eqs. (|A.1[) , ( |A.2| ) and (kA)), these potentials can be 
expanded near the black hole singularity in terms of the tortoise coordinate, respectively, as 



Vr(r) ~ -- 2 



Vy{r) ~ -— 2 



V S (r) 



1 

Tz 2 



1 o 2 w sl ' : 
i-2 2 -W S v(-^~ 



where 



WW 



1 



-2(d - 2)m]d r 2 
1 

-2(d - 2)m]^2 



-2(d - 2)m]a=* 



2(d-3) 2 41(1 + d - 3) 
2d -5 + d-2 

2(d 2 - 8d + 13) 4£(£ + d - 3) 



2d - 15 



+ 



d- 2 



2d 3 - 24d 2 + 94d - 116 4(d 2 - 7d + 14)[^ + d - 3) - d + 2] 



(2d-5)(d-2) 



+ 



(d- l)(d-2) 5 



and we have rescaled the tortoise coordinate z = ur* 
Near r = oo one has 



V T (r) ~ 0, 
WW ~ 0, 
Vs(r) ~ 0. 



The RN Case: The background space-time metric is 



t( \ 1 2m 



Near the black hole singularity (r ~ 0), the tortoise coordinate can be expanded as 

dr 1 r 2d ~ 5 2m r 3d ~ 8 



f(r) 2d -5 q 2 3d - 8 q A 

where we have kept the second term in the expansion of r and have chosen the integration constant so 
that r* = at r = 0, and we have assumed r/r^ <C 1. Using this formula, one may obtain an approximate 
expression of r in terms of the tortoise coordinate r* 



\d-3 



d-2 



^ ~ [(2d - 5)g 2 (r -) 5 - M r,]^ - [(2d - h)q 2 (r^ 2d r^. (A.6) 
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As one applies (|A.6|) to the five different potentials (eqs. ( |A.1| ), ( |A.2| ) and (1A.4[ )), these potentials can be 
expanded near the black hole singularity in terms of the tortoise coordinate, respectively, as 



Vr(r) 



Vy±{r) 



1 

Iz 2 

1 

Iz 2 
1 

'Jz 2 



1 - 3T - W RNT 



z 



(r ) u> 



d-3 
2d -5 



1 - Jy± ~ W RNV± 



z 



1 - j s ± - W RNS ± 



where 



d-3 



z 

3d -7 



d-3 
2d— 5 



d-3 
2d-5 



3V± 



2-jr, 



2d- 5' 2d -5 

and we have rescaled the tortoise coordinate z = wr„. The expressions of Wrnt, W rnv ± and Wjyy'^i in 
these equations are 



W. 



W RNT = 0, 

4m(d - 2) (d 2 - d - 4) - 4(3d - 8)[(d - l)m T A] 



(2d -5) (3d -8) 

W0ws± = -[2m(d - 2)(4d - 10)g 2 + W±](r ) d - 3 [(2d - 5)q 2 ] ^ , 



(r ) d - 3 [(2d-5)g 2 ] 



2-d 
2d-5 



where 



W 4 



-mg 2 [-(d - l)(3d - 8)(1 - J2) + 2(4d 2 - 15d + 12)] - 2(d - 2)(3d - 8)mg 2 
4(3d - 8) \l{t + d - 3) - (d - 2)]g 4 



(d-l)(l-n)m 
-mg 2 [(d - l)(3d - 8)(1 - 0) + 2(d - 2) 2 ] - 2(d - 2) (3d - 8)m<? 2 
4(3d - 8)[£(^ + d - 3) - (d - 2)]q A 



(d-l)(l + 0)m 
with the definition, here and below, that 



n = 

For the extremal case, q — > m, we have 

F T (r) - 

Vv±(r) ~ 

Vs±(r) ~ 



4[^ + d-3) -d + % 2 



1 + 



1 

"4^2 
1 

1 

4? 



(d — l) 2 m 2 



l — i 2 — w ex 
1 " " 



d-3 1 

Z \ 2d- 



2 , « -< r / "- \ 2d=5 



d-3 ~] 

Z \ 2d-5 
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where 



W e R X NT = 0, 

W e X 4m(d - 2)(ri 2 - d - 4) - 4(3ri - 8)[(d - l)m T A ex '] , M 

" (2d-5)(3d-8) [( " ' J 



and 



J{d - if (d -3) 2 + 2{d- 2) (d-3)(e(£ + d-3)-(d- 2) )m. 



Near r = oo one finds 



Vr(r) ~ 0, 
WW ~ 0, 
Vs(r) ~ 0. 



The Schwarzschild dS Case: In this case, the black hole is uncharged, so one has Vy+ = V\j- and Vg+ = Vg- 
For convenience, we denote them by Vy and Vs, respectively. The background space-time metric is 



f(r) = 1 



2m 



Ar 2 , 



with A > 0. Near the black hole singularity (r ~ 0), the tortoise coordinate can be expanded as 

(l r l j.d—2 -y r 2d-5 



f(r) d- 2 2m 2d - 5 (2m)' 



+ 



where we have kept the second term in the expansion of r and have chosen the integration constant so 
that r* = at r = 0. Using this formula, one may obtain an approximate expression of r in terms of the 
tortoise coordinate r* 

{d-2)u 



-2m(d - 2)r*]<*-2 + 



2d -5 



(A.7) 



As one applies ( |A.7| ) to the three different potentials (eqs. (|A.1[) , ( |A.2| ) and (A.4)), these potentials can be 
expanded near the black hole singularity in terms of the tortoise coordinate, respectively, as 



V T (r) 
Mr) 



1 

"4? 
1 

li 2 



1 - 2 - WsdST 



V s (r) ~ --L 



d-3 n 

2 \ d-2 



i -> 2 wwP l! 

\UJ. 
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where 



WsdST 

Wsdsv 
Wsdss 



[-2(d- 2)m}~ 
1 

[-2(d- 2)m]3=s 
1 



2(d-3) 2 41(1 + d - 3) 
2d - 5 + d- 2 

2(d 2 -8d+13) 4^(£ + d-3) 



2d - 15 



+ 



d- 2 



2d 3 - 24d 2 + 94d - 116 4(d 2 - 7d + 14)[£(^ + d - 3) - d + 2] 



2(d-2)m]^2 L (2d-5)(d-2) 



(d- l)(d- 2) 5 



and we have rescaled the tortoise coordinate z = cor*. It is easily seen that this is just like in the pure 
Schwarzschild case. 

Near r = oo one finds 



which leads to 



Vj(r) 
Vy(r) 
V S (r) 



1 f dr 1 
r*lr]~— x J ^ = r* + y r , 



d(d-2) cu[(d-l) 2 -l] 



4(r, 



a(z - zy 



(d-2)(d-4) _ w[(d-3) 2 -l] 



4(r* - r*) 2 4(z - z) 2 

(d-4)(d-6) cu[(d-5) 2 -l] 



4(r* 



4(z - z) 2 



as one defines z = ur* and z = wr, . 

The RN dS Case: The background space-time metric is 

2m 



f(r) 



+ 



yd— 3 y>2d—Q 



Ar 2 , 



with A > 0. Near the black hole singularity (r ~ 0), the tortoise coordinate can be expanded as 

dr 



/ 



I r 2d—5 1m r^d— 8 

+ » ; t ; + 



/(r) 2d -5 g 2 3d -8 g 4 



where we have kept the second term in the expansion of r and have chosen the integration constant so that 

r* = at r = 0, and we have assumed v/tq <C 1 ( = — -y/ m 2 — q 2 ^j d 3 represents the inner horizon 
of the RN black hole) . Using this formula, one may obtain an approximate expression of r in terms of the 
tortoise coordinate r* 



— ~[(2d-5)g 2 (r 



\5-2<f„ 



l^br 2m ( r 



\d-3 



(3d - 8)g 2 



[(2d-5)<7 2 (r 



\5-2d 



d~2 

n 1 2d ~ 5 ■ 



(A.8) 
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As one applies (|A.8|) to the five different potentials (eqs. ( |A.1| ), ( |A.2| ) and ( A.4 )), these potentials can be 
expanded near the black hole singularity in terms of the tortoise coordinate, respectively, as 



V T (r) 
Vy±(r) 
V 5 ±(r) 



1 

1 

"4^2 

1 

'Iz 2 



1 - 3t ~ ^RNdST 



[ r -)2d- 5uj 



d-3 
2d-. r > 



1 _ Jy± — W RNdS v± 



1 - is± ~~ W RNd $s± 



z 


(r ) 2d - 




z 




; r -)M- 





d-3 
L'd-r, 



d-3 

2d-5 



where 



d-3 . 3d -7 „ . 
Jt=3s±=tti — 7> Jy± = — — ^ = 2-jt, 



2d -5 2d -5 

and we have rescaled the tortoise coordinate z = wr». The expression of WRNdST > Wijjvdsv* an d WijiVdSSi 
in these equations are 



ir, 



WWdST — 0, 

4rn(d - 2)(d 2 - d - 4) - 4(3d - 8)[(d - l)m T A] 



RNdsv± 



\d-3 



(2d-5)(3d-8) v 
W RNdSS ± = -[2m(d - 2)(4d - 10)g 2 + W ± ](r [(2d - 5)g 2 j 



o > 

7-3d 



[(2d-5)g 2 



where 



W 4 



-mg 2 [-(d - l)(3d - 8)(1 - fi) + 2(4d 2 - 15d + 12)] - 2(d - 2) (3d - 8)mg 2 
4(3d - 8) [£(£ + d - 3) - (d - 2)]q A 



(d-l)(l-n)m 

-mg 2 [(d - l)(3d - 8)(1 - fi) + 2(d - 2) 2 ] - 2(d - 2)(3d - 8)mq 2 

4(3d - 8)[g(l + d - 3) - (d - 2)}g 4 
(d- l)(l + fi)m ' 



This is just like the pure RN case. 
Near r = oo one finds 



1 f dr 



1 

A? 



which leads to 



Vj(r) 
W(r) 
^s(r) 



d(d-2) w[(d-l) 2 -l] 



4(r* - r*) 2 4(z - z) 2 ' 
(d-2)(d-4) _ w[(d-3) 2 -l] 



4(7%,- r*) 2 4(z-z) 2 
(d-4)(d-6) w[(d-5) 2 -l] 



4(r„ - r*) 2 4(z-z) 2 
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as one defines z = ur* and z = wr t . 



The Schwarzschild AdS Case: In this case, the black hole is uncharged, so one has Vy+ = Vy- and Vg+ = Vg- 
For convenience, we denote them by Vy and Vs, respectively. The background space-time metric is 



f(r) 



1 



2m 

r d—3 



Ar 2 , 



with A < 0. Near the black hole singularity (r ~ 0), the tortoise coordinate can be expanded as 

(l r l j.d—2 -y r 2d-5 



d - 2 2m 2d - 5 (2m) 2 



+ 



where we have kept the second term in the expansion of r and have chosen the integration constant so 
that r* = at r = 0. Using this formula, one may obtain an approximate expression of r in terms of the 
tortoise coordinate r* 



r ~ [-2m(d - 2)r*]3^ + (d 2)r * , 

la — 5 



(A.9) 



As one applies ( |A,9| ) to the three different potentials (eqs. ( |A.l ), ( |A.2| ) and (A.4)), these potentials can be 
expanded near the black hole singularity in terms of the tortoise coordinate, respectively, as 



V T (r) 
Vy(r) 
V s (r) 



1 

Tz 2 
l 

Az 2 
1 

Tz 1 



d-3 

1 " 2 - WsAdST (- '' J 

1-2* -W SA dsvP H 

l-tf-WsAdSst-^ 1 

\L0 



where 



WsAdST 
WsAdSV 
WsAdSS 



1 



-2(d - 2)m]d=2 
1 

-2(d - 2)m]^2 
1 

-2(d - 2)m]^2 



2(rf-3) 2 41(1 + d - 3) 
2d - 5 + d-2 

2(d 2 -8d + 13) 4£(£ + d-S) 



2d -15 



+ 



d-2 



2d 3 - 24d 2 + 94d - 116 4((f 2 - 7d + 14) [i(* + d - 3) - d + 2] 



(2d-5)(d-2) 



(d-l)(d-2) 2 



and we have rescaled the tortoise coordinate z = wr t . It is easily seen that this is just like in the pure 
Schwarzschild case. 

Near r = oo one finds 



1 f dr 



1 

A I r' 
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which leads to 



Vj(r) 
W(r) 
V S (r) 



d{d-2) w[(d-l) 2 -l] 



4(r* - r*) 2 4(z - z) 2 

(d-2)(d-4) _ u;[(d-3) 2 - 1] 
4(r* - f*) 2 4(z - z) 2 

(d-4)(d-6) _ u;[(d-5) 2 -l] 
4(r*-f*) 2 



as one defines z = cur* and 2 = wr„ . 

The RN AdS Case: The background space-time metric is 



f(r) = 1 



2m 



+ 



4(z - z) 2 



Ar 2 



with A < 0. Near the black hole singularity (r ~ 0), the tortoise coordinate can be expanded as 

+ ^ r :— + ••■ • 



/(r) 2d -5 g 2 3d -8 q 4 
where we have kept the second term in the expansion of r and have chosen the integration constant so that 

r* = at r = 0, and we have assumed r/r^ <1 ( r$ = im— \J m 2 — q 2 



i 

d-3 



represents the inner horizon 
of the RN black hole) . Using this formula, one may obtain an approximate expression of r in terms of the 
tortoise coordinate r* 



[(2d-5)(? 2 (r - 



\5-2d 



2m (r 



\d-3 



■[(2d-5)^(r -) 5 - M r*]^. 



(A.10) 



(3d-8)g 2 

As one applies ( A,10| ) to the five different potentials (eqs. (|A.1[) , ( |A.2| ) and ( |A.4j) ), these potentials can be 
expanded near the black hole singularity in terms of the tortoise coordinate, respectively, as 



V T {r) 
Vy±{r) 
V s ±(r) 



1 

'Iz 2 
1 

'Iz 2 

1 

'Iz 2 



1 — 3t ~ WrnAcLST 



z 



(r ) 2d 5 w 



d-3 ■ 
2d-5 



1 - 3v± - W RNAdS y± 



z 



\2d-5 



where 



3T = J5± 



1 - 3s± ~ W RN AdSS± 
d-3 



{r^) 2d -^u 
3d -7 



d-3 ■ 
2d-5 



d-3 

2d-5 



2-iT, 



2d -5 K 2d -5 

and we have rescaled the tortoise coordinate z = ujr*. The expression of W R NAdST-, ^RNAdSV ± an d 
W RN AdSS ± m these equations are 

WRNAdST = 0, 



W, 



RNAdSV ± 



Am(d - 2)(d 2 - d - 4) - 4(3d - 8)[(d - l)m T A] 
(2d-5)(3d-8) ( 

= ~[2m(d - 2)(4d - 10)g 2 + W±](ro )^ 3 [(2d - 5)g 2 



2-d 

\d— 3 IVoJ K\„21 2d-5 



[(2d-% 2 
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where 



W+ = -mq 2 [-(d - l)(3d - 8)(1 - fi) + 2(4d 2 - 15d + 12)] - 2(d - 2) (3d - 8)mq 2 - 

4(3d - 8) [1(1 + d - 3) - (d - 2)]g 4 
(d- 1)(1 - fi)m ' 

W_ = -mq 2 [(d - l)(3d - 8)(1 - fi) + 2(d - 2) 2 ] - 2(d - 2)(3d - 8)mg 2 - 

4(3d - 8) [1(1 + d - 3) - (d - 2)]g 4 
(d- l)(l + fi)m ' 



This is just like the pure RN case. 
Near r = oo one finds 




which leads to 



V T (r) 



d(d - 2) 



W [(d-1) 2 -1] 



W(r) 



V s (r) 



4(r* - r*) 2 

(d-2)(d-4) 
4(r* - f*) 2 

(d-4)(d- 6) 
4(r* — f*) 2 



4(z-z) 2 



w[(d-3) 2 -l] 



cu[(d-5) 2 -l] 



4(z - zf 



4(z - z) 2 



as one defines z = wr» and z = tor* . 
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